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Abstract 



The post Gaussian effective potential m D = 2+2e dimensions is evaluated 
for the Ginzburg-Landau theory of superconductivity. Two and three loop 
integrals for the post Gaussian correction terms in D = 2 + 2e dimensions are 
calculated and e-expansion for these integrals are constructed. In L) = 2 + 2e 
fractal dimensions Ginzburg Landau parameter turned out to be sensitive 
to e and the contribution of the post Gaussian term is larger than that for 
D = 3. Adjusting e to the recent experimental data on n(T) for high -Tc 
cuprate superconductor Tl2Ca2Ba2CusOio{T£ — 2223), we found that e = 
0.21 is the best choice for this material. The result clearly shows that, in 
order to understand high - Tc superconductivity, it is necessary to include the 
fluctuation contribution as well as the contribution from the dimensionality 
of the sample. The method gives a theoretical tool to estimate the effective 
dimensionality of the samples. 
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I. INTRODUCTION 



The Ginzburg-Landau (GL) theory of superconductivity [1] had been proposed long 
before the famous BCS microscopic theory of superconductivity was discovered. A few years 
after the appearance of the BCS theory, Gorkov derived the GL theory from the BCS theory 
[2]. Since then, the GL theory has remained as a main theoretical model in understanding 
superconductivity. It is highly relevant for the description of both type -I [3] and type II 
superconductors, even though the original BCS theory is inadequate to treat both materials. 
The success of the GL theory in the study of modern problems of superconductivity lies on its 
universal effective character in which the details of the microscopic model are unimportant. 

Even at the level of meanfield approximation (MFA), the GL theory yields significant 
information such as the penetration depth (£) and the coherence length (^) of the supercon- 
ducting samples. Many unconventional properties of superconductivity connected with the 
break down of the simple MFA has been studied both analytically [4] and numerically using 
the GL theory [5] . Particularly, the fluctuations of the gauge field were studied recently by 
Camarda et. al. [6] and Abreu et. al. [7] in the Gaussian approximation of the field theory. 
The effective mass parameters of the Gaussian effective potential (GEP), Q and A , were 
interpreted as inverses of the coherent length ^ — 1/Q and of the penetration depth i — 1/A, 
respectively. 

In our previous paper [8] we have estimated corrections to the Gaussian effective potential 
for the U{1) scalar electrodynamics, which represents the standard static GL model of 
superconductivity. Although it has been shown that the correction is significant in D — 3 
dimensions, it was not large enough to explain the experimental findings. At the same time, 
we have investigated the role of quasi two dimensionality in the high Tc superconductivity, 
by calculating the Gaussian effective potential for D = 2 + 2e. It was found that the 
dimensional contribution at the Gaussian approximation level gives the correction in the 
right direction, but is not large enough to explain the experimental data [8]. However, it 
is known that fluctuation contributions are much larger in lower dimensions. Therefore, 
it is necessary to investigate whether the post Gaussian correction terms in D = 2 + 2e 
dimensions provide significant contribution to the mean field result, in order to understand 
the layered structure of the high Tc superconductivity. In the present paper, we study the 
role of the post Gaussian contributions in D = 2 + 2s dimensions by using the method 
developed in [8]. 

The paper is organized as follows: in Section II the GL action is introduced and basic 
equations are derived; in Section III, the theoretical results for D = 2 + 2e will be compared 
to existing high Tc experimental data, so that the role of fractal dimensions can be discussed. 
In the Appendix we calculated two and three loop integrals in T* = 2 + 2^ dimensions. 
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II. BASIC EQUATIONS FOR THE EFFECTIVE MASSES 



The Hamiltonian of the model and cxphcit expressions for the effective potential in 
Euclidean D-dimensional space were given in [6-8]. Here we bring the main points for 
convenience. The effective potential, i.e., the free energy density, V^q = T jV is defined as 

l-eff=-lnZ (2.1) 

where the partition function is 

Z = ^ V(j)VATexp{- J d^xH + J d^xjcj) + (jaA)}. (2.2) 

The Hamiltonian density is given by 

= i(V X Ar + + \m^' + + \e^''A' + ^(VA)^ (2.3) 

where we have introduced a gauge fixing term with the limit — being taken after the 
calculations are carried out. Note that, we are using natural units employing (coherence 
length at zero temperature) and (critical temperature) as the length and the energy scales, 
respectively, introduced by [9]: 

° (2.4) 

Using the method introduced in refs. [8,10,11] one finds following effective potential: 

V^Q^VG + i^Va (2.5) 

where Vq is the Gaussian part: 

Vg = hiSl) + i/i(A) + ]^mm + Xct>l + ho{.^)[m' + 6XIo{n) + UX^] (2.6) 
+ 7o(A)[-A^ + e%(Q) + e20g], 

and AVg is the correction part: 

^Vg = [-^e%(A) - 18l2(l^)A2]0o^ + {-3Xh{n)[-n^ + m^ + 2Io(A)e2 + 12AIo(fi)] 

- e%(A)[-A2 + e^Io(Q)] - 8A%(Q, Q) - ^e%(A, Q)}(/.o' - + m^ + 2Io(A)e2 

o o 

+ 12AIo(n)]2 - ll2(A)[-A2 + e%(0)]2 - ^e%(A,n) - \x\^{^,^). (2.7) 
In the above following integrals are introduced: 
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2 f d^k 
''^ ' ~ (27r)-D J (P + M2)2' 



/3(Mi,M2) 

h{Mi,M2) 



^_ I d^kd^P (2.8) 



(27r)2^ J (A;2 + M2)(p2 + Ml){{k + pf + Mf^ ' 
d'^kd^pd^q 



(27r)3^ 7 (A;2 + M2)(p2 + M2)(g2 + Ml) {{k+p + qf + Mf) ' 

For D = 3 — 2e, these integrals were calculated in dimensional regularization in ref. [12] and 
for D = 2 + 2e in the Appendix of the present paper. 

The parameters Vt and A are determined by the principle of minimal sensitivity (PMS): 

9KfF A2 A "2 



- 0.0362— + -[0.075(J^2 - m^) - 0.108A(ln ^ + 1) - 0.911(/)§A 



' £2 ' • V ' ■ y 

2 2 

+ 0.145A2 In^ ^ + [(0.290 + 1.823(^g)A2 + 0.151A(m2 - O^)] In ^ + 0.064A2(1 + 3^^)^ 

+ {n? - J^2)[0.039(m2 - f^^) + 0.954A^^0.151A] - O.IOSA^ In^ ^ 

„2 _ 

+ [0.113(^7^ _ m'^)\ - a2(0.326 + 1.367(/>^)] In^ !t- + m'^ - m^)iQ.9bAX(Pl (2.9) 
+ 0.227m2)A - 0.039(0^ - rr?)) - A2(0.133 + 5.729(/)^ + 3.215(^§)]ln ^ 

- A^(O.96O0g + 5.72^^ + 0.144) + £[{^'^ - m2)(0.954A^g - 0.062A - 0.039(O^ - rr?)] + 0{e^) = ; 
^ = (0-334O2-0.319A) ^^^^ _ q 334^2 ^ 3^9;^ _ o.33402)ln + (4.015A - 4.205^)^)^2 

,,2 

- 1.0031)2 + 0.334m2 + £{(0.1671)^ - 0.479A) In^ 

i)2 

2 2 

+ [(0.3341)2 _ o.639A)ln ^ + 1.0031)2 - 0.334m2 - 4.015A^g]ln ^ 

/,2 

+ [(4.2051)2 - 4.O15A)0^ + 0.3341)2 _ o.334m2]ln ^ - 0.262A - 4.9431)2 A2 (2.10) 
+ 8.4101)2,^2 ^ 0.2751)2} + 0(£2) = 0, 

where we denote optimal values of 1) and A by and A, respectively, and 0o is a stationary 
point defined from the equation: 

agff _ AO^A^(0.450A-0.477fi^) _^ ^^^^^^^ _ ^ _ ^^^^^^ _ 

+ 0.477(1)2 - m2) - 0.119A + — ^ + {[0.683A - 0.2381)2] ln2 ^ (2.11) 

„2 

+ [5.729A0^ + 0.477(m2 - 1)2) + 0.239A]ln ^ + 0.240A - 0.3921)2}e + 0(£2) ^ g. 
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In the equations (2.9) - (2.11) we have used e expansion of the loop integrals exphcitly and 
numerical values of ^o, Tc and e. For the cuprate Tl2Ca2Ba2Cu30io{Ti — 2223) these values 
are 

Co = IMnm, Tc = 121.5E:, = IQ-KakBTcio/fic = 0.0000264. (2.12) 

III. RESULTS AND DISCUSSIONS 

The solutions of the Eqs. (2.9) - (2.11) are related to the experimentally measured GL 
parameter k as k — ij^ — ft/ A. We make an attempt to reproduce recent experimental 
data on k,{T) [13] for high -Tc cuprate superconductor Tl2Ca2Ba2Cu30io{Ti — 2223). 

For this purpose, we adopt usual linear T dependence of parametrization of m and A as: 

= ml{l - t) + Tml, A = Ao(l - r) + tAc, t = T/Tc, (3.1) 

and calculate k by solving nonlinear equations (2.9) - (2.11). Due to the parametrization 
(3.1), the model has in general six input parameters: mg, Aq, m^, Ac, (coherent length) 
and Tc (the critical temperature). The experimental values for the cuprate Tl — 2223 are 
Co = 1.36nm and Tc = 121. SX. To determine other four parameters we used the following 
strategy. For each given £, the parameters ttIq and Aq are fitted to the experimental values 
of C and £ at zero temperature: Co — 1.36nm, — 163nm. In dimensionless units, (2.4), 
we have Qq = ^(t = 0) = 1 and Aq = A(t — ^) — ^o/h — 0.0083 which are used to 
calculate and Aq from the coupled equations (2.9) - (2.11). This procedure gives the e 
dependence of rriQ which is presented in Fig. 1 (sohd hne) . As in the case of the Gaussian 
approximation [8], remains positive only for very small values of £, although nonlinearity 
produces several rm? — Q solutions in this case. We believe that this smallness again indicates 
the reliability of the present post Gaussian approximation method. 

The parameters ml and Ac are fixed in the similar way for each given e. Actually the 
quantum fluctuations shift ml from its zero value given by MFA. On the other hand, the 
exact experimental values of m^ and Ac are unknown, since the GL parameter at T = Tc is 
poorly determined. For this reason, we used the experimental values of Cc and ^c at very close 
points to the critical temperature, Tc = 0.98 which corresponds to VL^ = CI{tc) = 1/Cc = 0.128 
and Ac = A(rc) = 1/4 = 0.0043 (kc = 29.6). Then solving the equations (2.9) - (2.11) 
numerically with respect to mc and Ac, wc fix these parameters. 

After having fixed the input parameters, the temperature dependence of CI{t), A(r) 
as well as the GL parameter k = Q{t)/A{t) arc established by solving the gap equations 
numerically for each e. Clearly, the solutions of nonlinear gap equations are not unique. In 
numerical calculations we separated the physical solutions by observing the sign of 0q, which 
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should be positive and that the effective potential at the stationary point V^gp(0o) should 
has a real minimum at this point. For e > 0.1, there is a possibility to adjust e to the recent 

experimental data on k,{T) [13] for high -Tc cuprate superconductor Tl2Ca2Ba'2CuzOn){TC — 
2223). Our calculations show that, the best choice of e is found to be e = 0.21. The 
appropriate k{t) is presented in Fig. 2 (solid line). The dashed line in this figure shows k{t) 
for D = 2>. This fitting process allows us to get an estimation on the effective dimensionality 
of the high - superconducting materials. 

IV. SUMMARY 

In the present paper, we have carried out two and three loop calculations on the Ginzburg- 
Landau effective potential beyond the Gaussian approximation for D = 2 + 2£: fractal di- 
mensions. The result clearly shows that the higher order corrections are substantially large 
to explain the existing experimental data. 

This result strongly suggests that in order to explain the experimental data on high - 
Tc superconductivity it is necessary to include the fluctuation contribution as well as the 
contribution from the quasi two dimensionality. We have found that the GL parameter is 
rather sensitive to e when the loop corrections to the simple Gaussian approximation are 
taken into account. The optimal value of e for the cuprate {Tl — 2223) is £ = 0.21. It would 
be interesting to estimate optimal e in fractal dimensions for other cuprates also. 

It is to be noted that we have calculated two and three loop integrals m. D — 2 -\-2e 
dimensions using the method of dimensional regularization. 
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Appendix 

EXPLICIT EXPRESSION FOR THE LOOP INTEGRALS IN D = 2 + 2£ 

DIMENSION. 

Here, we consider the loop integrals defined in Eqs. (2.8) m. D = 2 + 2e dimensions. In 
dimensional regularization the integrals loiiri), Ii{m) and hiin) can be easily calculated in 
momentum space: 
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(27r)^ (m2 + p2) ^ ' r(L»/2)(27r)^ Jo (A;^ + m^) 
with X = //^/m^. 

Two and three loop integrals {I3 and 74) require a little more effort. It is more convenient 
to evaluate them in coordinate space rather than in momentum space, since 

1 r dPkd^pd^q 1 

74(^1, M2) - j—^ J ^ 2 ^ 2 ^ 2 ((k + v + 0)2 + Ml) ^ ^ ^ 



(27r)3^ y + Mf) (p2 + Mf) (g2 + m|) {{k+p + q)^ + Ml) 
'^rJd^rGl{r)Gl{r), 



where G'„(r) is the Fourier transform of the propagator l/(k + M^) (n — 1, 2): 

^^'^ = I (2.)^(P + m2) = (A.3) 

and K^,[z) is the modified Bessel function. In dimensional regularization, for D = 2 + 2£, 
G{r) is simplified as 

G{r) = {^r^^^Umr). (A.4) 

Now, substituting (A.4) into (A.3) one notices that unlike in the case of = 3 — 2^ , in 
D = 2 + 2e dimensions there is no singularity at small r and hence the integration can be 
performed directly from r = to r = 00 without splitting radial integration into two regions 
with small r and large r. 

The case with equal masses. Mi = M2 = m, can be done analytically: 



'^^"^^ = (2.)^-im2r(i + ,) -^^(^)' ^^(^) = i t'^'-r-Ksrat 



(A.5) 



for = 3, 4, where the integrals /3(£) and h{e) are expressed in term of the hypergeometric 
functions: 



7 



/3(^) { r(3/2-£) ^^^f^' 1 - 2 " 4^ 



-2r(l-£)2Fi[l,l-£;|^]}, (£<0.5); 



~ , , r(£)r(l -£) ,eT^{-e) ^, 13 



4^(2£-l) ^ ^2' ' '2 ' ' ^ 



4-V^r(i-35)r(i + £)r(-2£) ^ i ^3 



(A.6) 



r(i-2e) 2F^[l-3e,--e■,--2e■,l]}, < 1/3). 

The method of ref. [14] gives the following e expansion: 

Um) = ^ I J 0.5917 + £(0.6629 - 1.1835 In x) + Ofe^)], 
IJm) = ^ J l.188 - e(2.759 + 3.5656 In x) + Ofe^)] 

which is used in our practical calculations. 

The case with nonequal masses is rather complicated and cannot be done analitically 
in general. However, in the particular case, when a = M2/M1 < 1 ^ the problem may 
be overcome by expansion in power series in a. We shall illustrate this approximation for 
Is{Mi,M2) below. Using Eq.s (A. 3) and (A.4) one obtains 

where 

_ roo 

I^{a,s)= tKHt) {at)-'K,{at). (A.9) 
Jo 

Now using the series expansion of ^",^(2;) 

^^^'^^ 2 ^W^) W^) ^ 



in the present paper a = 1/k where k 80 in the large range of temperature 
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one may expand the factor (at) ^K^{at) in power series of a and integrate (A. 9) analytically 
to obtain: 

- 3q;-'"[4£ - 6 + q;'(2£ - 1)] + 0{a^)}. (A.ll) 
Inserting Eq. (A.ll) into the Eq. (A. 8) one obtains the following e expansion : 

/3(Mi,M2) = ^ {108(l-lna) -3a^(61na-5) +£[a^(-181n^a + (361na;i + 18)lna 
oD47r ikZ]^ 

- 30 In xi + 4) - 216 In xi - 108 In^ a + 216 + 216 In xi In a + 0{e^)]}. (A.12) 
Similarly, one may calculate hlMi, M2) to obtain it's e expansion: 

/4(Mi,M2) = o , ^2 {4a^(2 + 9 In^ a - 6 In a) - 108 In^ a + 190.9588 In a - 280.5109 

172877'^ 

+ £ [q;^(72 In^ a - (60 + 108 In xi) In^ a + (72 In xi - 28) In a - 24 In xi + 23.4519) 

- 3601n^ a + (547.8351 + 3241nxi)ln2 a - (572.87641nxi + 337.6413) Ina 

+ 841.5330 Inxi - 806.1519 + 0{e^)]} (A.13) 

where, for simplicity, we used explicit values of constants such as 7, C(3), ln(2), etc. 
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FIG. 1. The parameter of the GL model v.s. e in fractal dimension = 2 + 2e. The solid 
and dashed lines are for the temperatures T = (|3ind T = O.GTc respectively. 



Post Gaussian D=2+2e 

— Post Gaussian D=3 
■ Experiment 




FIG. 2. The GL parameter, k, in D = 2 + 2e (solid line) and in D 
calculated in the PostGaussian approximation. 12 



= 3 (dashed line) cases 



